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Abstract. The method of (/-oscillator lattices, proposed recently in hep-th/0509181 , pro- 
vides the tool for a construction of various integrable models of quantum mechanics in 2 + 1 
dimensional space-time. In contrast to any one dimensional quantum chain, its two dimensional 
generalizations - quantum lattices - admit different geometrical structures. In this paper we 

o 

, consider the g-oscillator model on a special lattice. The model may be interpreted as a two- 

, dimensional Bose gas. The most remarkable feature of the model is that it allows the coordinate 

Bethe Ansatz: the p-particles' wave function is the sum of plane waves. Consistency conditions 
is the set of 2p equations for p one-particle wave vectors. These "Bethe Ansatz" equations are 

(N ■ 

. the main result of this paper. 



d 



1. Introduction 



Integrable models of quantum mechanics in discrete space-time describe a system of interact- 
ing quantum observables (spin operators, oscillators, etc.) situated in sites of an one dimensional 
chain (the models in 1 + 1 dimensional space-time), or in the vertices of a two dimensional lattice 
(models in 2 + 1 dimensional space-time). The subject of this paper is the second variant - the 
quantum integrable lattices. 

Formulation of a realistic integrable model of quantum mechanics in 2 + 1 dimensional space- 
time was a long-standing problem in the theory of integrable systems. Examples of higher 
dimensional integrable systems are known, but the integrability seemed to be a very high price 
j> \ paid to the detriment of a physical interpretation. 

•i-H . 

^ ■ Both known classes of integrable quantum lattices may be specified by the definition of the 

algebra of observables. In contrast to the quantum chains, where the algebra of observables may 
be e.g. an evaluation representation of any quantum group, in the case of quantum lattices the 
algebra of observables may be either the local 1 Weyl algebra or the local g-oscillator algebra 
0. Other examples of integrable quantum lattices are not known. 

When the algebra of observables (and its representation) is fixed, there exists an almost 
unique way to produce the key notion of the integrability - a complete set of commutative 
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quantum chains case 



^The term "local" means that the quantum observables in different vertices commute, as well as it is in the 
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operators. For the spin chains it is the way of Lax operators and transfer matrices with e.g. 
fundamental representations in the auxiliary space. In both classes of quantum lattices the 
integrals of motion may be produced by a decomposition of a determinant of a certain big 
operator- valued matrix with respect to two spectral parameters 013] • Therefore the integrals of 
motion rather have a combinatorial nature. Both classes of quantum lattices allow a definition of 
another kind of operators: the layer-to-layer transfer matrices, related to intertwining relations, 
tetrahedron equations etc. Such representation-dependent layer-to-layer transfer matrices are 
the kind of discrete-time evolution operators for the quantum lattices, their structure is much 
more complicated, but since their integrals of motion are known, it is not necessary to consider 
them at the first step. 

The concept of quantum lattices implies a certain feature which is not applicable to the 
quantum chains. One may vary not only a size of the lattice, but its geometry and topology as 
well. This is the reason why we discuss the classes of models. The methods of Weyl algebra 
and g-oscillator algebra are the frameworks for construction of various quantum lattices. 

Mainly the simplest geometry was considered before - the square lattices with periodical 
boundary conditions 2 . The simple square lattice with the sizes n x m may be regarded as the 
length- m chain of its length-n lines, and therefore such quantum lattice is effectively a quantum 
chain. In particular, the quantum lattice with the Weyl algebra at a root of unity corresponds 
to auxiliary transfer matrices 1] of the U q {gl n ) generalized chiral Potts model [HI El- The q- 
oscillator quantum lattice corresponds to the auxiliary transfer matrices for reducible oscillator 
representation of U q (gl n ) 2 . From the point of view of quantum mechanics in 2 + 1 dimensional 
space-time, the square lattices provide a non-realistic models since the translation operators do 
not belong to the set of integrals of motion and therefore momenta of eigenstates can not be 
defined. In the quantum chain interpretation, it means that U(l) charges are the variables, and 
the chain is homogeneous only on subspaces where all U(l) charges are the same in all sites. 
Otherwise, the translation invariance is lost. 

The issue is evident. One has to consider a lattice which is not a chain of its lines. We define 
such lattice in the next sections. The framework of g-oscillator model on our special lattice 
immediately produces a kind of two dimensional Bose gas with physical dynamics. The Fock 
vacuum is the natural reference state, the total ^-oscillators occupation number - the number 
of bosons - is conserved, and the p-particles wave function is the superposition of one-particle 

2 Examples of Weyl algebra non-square lattices were mentioned in 0. In particular, the Weyl algebra framework 
on the different shapes of the lattices produces the relativistic Toda chain as well as the quantum discrete Liouville 
model. 
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plane waves. Each plane wave is characterized by two components of its momentum, and the 
consistency conditions give 2p equations for p two-components momenta. These consistency 
conditions look much more complicated than the Bethe Ansatz equations for a quantum chain, 
we do not investigate them in details here. The aim of this paper is just to present the method, 
the model, and the "Bethe Ansatz" equations. 



2. Framework of (/-oscillator algebra. 
We start with the formulation of a generic ^-oscillator lattice. 

Let £ be a lattice formed by a number of directed lines on a torus. Pairwise intersections of 
the lines are the vertices of the lattice. Let the vertices are enumerated in some way, for instance 
by the integer index j = 1, 2, 3, . . . , A, where A is the number of vertices. 

The g-oscillator generators Xj,yj and hj are assigned to vertex j. The oscillators for the 
different vertices commute. Locally, we define the g-oscillator algebra by 



(1) 



1 _ (,2+211, 



1 



,2h, 



Consider next the system of non-self-intersecting paths along the edges of the lattice. The lines 
of the lattice are directed, and we demand that paths must follow the orientation of the edges. 
The paths may go through a single vertex in one of six variants shown in Fig^ The vertex 



_[ 



r 



Figure 1. The "weig hts" fj of j th vertex, j = 1, A. 

"weight" fj is associated with each variant of bypassing as it is shown in Fig. ^ Coefficients Xj, 
[ij and Uj are extra vertex C- valued parameters, they are related by = —q~ 1 \j[ij. In general, 
they are free up to a single condition: if two vertices j and j' are formed by the intersections of 
the same lines, their C-valued parameter must be the same jl]. 
For any path P let 

(2) t P = n f 3 . 

along P 

Note, if a path does not touch a vertex, then it vertex's contribution to tp is just the unity 
according to the leftmost variant of Fig. ^ Recall in addition, the ^-oscillator algebra is local, 
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elements of different vertices commute. Therefore the notion of the product along the path in 
© is well defined. 

Let A and B be the two basic homotopy cycles of the torus. Any path has a certain homotopy 
class, P ~ nA + mB. Let 

(3) t„ im = ^2 *>p . 

P : P~nA+mB 

Here the sum is taken over all possible paths of homotopy class nA + mB. Formula ends 
the formulation of (/-oscillator lattice's prescription. The point is that for any lattice C, formed 
by directed lines on the torus, the operators t n m form the commutative set: 

(4) ^n,m tn',m' — ^n',m' ^"n,m V n,m : n ,m . 

Moreover, if there are no lines of trivial homotopy class in C, the set of integrals of motion is 
complete. 

Integrals of motion may be combined into a polynomial of two variables, 

(5) T(«,t;) = ^u^ m t n , m . 

n,m 

Operator T has the structure of the layer-to-layer transfer matrix, equations (J3J are equiv- 
alent to the commutativity of transfer matrices with different spectral parameters u, v. The 
commutativity may be proven for simple lattices with the help of tetrahedron equation [HHHE]. 

3. The lattice 

Turn now to the definition of our special lattice. It is formed by the intersections of only two 
lines a and (5. With respect to the basic cycles A and B of the torus, the lines have the classes 

(6) a ~ NA - B , 13 ~ MB - A . 

Example of such lattice is given in Fig. [21 There the opposite dashed borders are identified 
(dashed lines are the cuts of the torus). The line a is directed up (along the ^4-cycle), the line 
(3 is directed to the right (along the B-cycle). The lines are intersecting in 

(7) A = NM-1 
points. 

We enumerate the vertices by number j, j G Za- The numeration is successive along the 
reverse direction of a-line, see Fig. |2 The shift j — ► j — 1 corresponds to a one-step translation 
along the a-line, whereas the sift j — » j + M corresponds to the one-step translation along the 
/3-line. 



ANSATZ OF HANS BETHE FOR A TWO-DIMENSIONAL BOSE GAS 



5 



T 

A 




B 



Figure 2. The lattice formed by the lines a ~ 4A - B and f3 ~ 3B - A. 



Now the lattice is formulated, and combinatorial rules of Fig. ^ ma y be applied. Note, since 
all the vertices of the lattice are formed by the same lines, all their C- valued parameters are the 
same, and one may put 



(8) 



A = fi = 1 , v 2 = — q 1 



Generating function © has the structure 



N-l M-l 



(9) 



T(n, v) = ^v- 1 ^ + v M u-^ + E E > 



n=0 m=0 



where TV is the total occupation number 



(10) 



Elements in Q are the values of both tjv-i and t_j Mj combinatorially they are the paths 
along all a or f3 lines (the third and the fourth variants of Fig. ^ f° r au vertices) . Elements to,o 
and tjv-i M-l correspond to the empty path and the complete path (the first and the second 
variants of Fig. [I\for all vertices) correspondingly, their values are 



(11) 



0,0 — 1 i tjV-l,M-l 



(-s-r . 
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In what follows, we will consider the elements to,i and t^o given by 
ti,o = -q- 1 52x d y j +M<? i+1+h * +a+ - +h * +1 '- 1 , 

(12) 

t 0jl = -q- 1 J2x j y j + 1 q h 3 + M +h 3 + 2M + - +h 3 + (N-l)M . 

3 

Combinatorial summand of to i is shown in Fig. |3 




B 



Figure 3. Path of the class B - an element of to i. 



4. ElGENSTATES. 

The total occupation number is the conserving quantity, and therefore we may construct the 
eigenstates of the model step-by-step just increasing the number of bosons for the Fock space 
representation of the g-oscillator algebra. The Fock vacuum |0) is defined by 

(13) Xj |0> = , hj |0> = VjGZ A . 

The Fock vacuum is evidently the eigenstate of the model, eigenvalues of all t n , m are zeros 
(except t 0l o, tjv-l,M-i and q^). 

In what follows, we concentrate on the diagonalization of to,i- Its eigenstates are defined 
uniquely. They will be evidently the eigenstates of t^o and, due to the uniqueness, they must 
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be the eigenstates of all the other integrals of motion. We will use below the normalized form 
of (HU): 

(14) r a = -^3^2*0,1* Tp = -JZ^* 1 ' ' 

4.1. One boson state. Consider the states with the single boson, TV = 1. Let 

(15) |0;> = Yj |0> 



be the basis of one-boson states. On this subspace both operators (fT4JI are the translation 
operators: 

(16) T a \(f>j) = \(f>j+l) , Tp \<f>j) = \4>j+M) ■ 

Therefore 

(17) |*> = \<t>i) , ^ = 1 , 

3 

is the eigenstate of the model, the eigenvalues of (|T1)) are 

(18) T a = |*) UJ- 1 , T p = |*) LU~ M . 

Since lo: lo a = 1 may take A different values, the set of (fT7|) is complete on the subspace J\f = 1. 

Formally, the spectrum of both translation operators Q18JI is defined by one parameter u. 
Note, eigenvalues a; -1 and u~ M are dual in the sense (uj~ m ) n = <j —1 . When N,M — > oo, one 
may talk about two independent components of momentum. Namely, let 

(19) k = N k a + k{3 , 
where fc a and /c^ are relatively small. Then, when N,M — > oo, 

^20) w _1 = g27r«fc/A _^ e 2nik a /M uj~ M = Q^iMk/A _^ Q 2nikp/N 

i.e. in the thermodynamical limit 2irk a /M and 2jrkp/N play the roles of independent compo- 
nents of the momentum. 

4.2. Two bosons state. Turn next to TV = 2. Define 

(21) Cj^ = Sj^+M + Sj t k+2M + 1" ^j,Jfc+(JV-l)Af • 

Then 

^ r a • yjVfc|0) = (q Ck '>y j+ iy k + q c i' k yjy k+1 )\0) , j^k, 

r a • y||0) = (l + g 2 )y fe+ iyifc|0) . 
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Consider a two-particles state with zero momentum of the mass center: 



(23) 



I*) = E ( *o(i + <z 2 rVf |o) + £ ^y^io)) 

j V fc>0 / 



Eigenvalue equation r Q |\&) = |\1/) r a reads in components 

r a ^ = (l + g 2 )*i , 

(24) 

in all other cases : T a ^k = ^>k-i + ^k+i • 
Periodicity condition is simply 

(25) V k = ^A_ fe . 
Define the double index: j = Mn + m — > j = (n, m), 

(26) ^Afn+m = *(n,m) , < m < M , 0<n<N . 

Equations (|24|) have in the general position the form r a ^/( n m ) = ^(^m-i) + ^( n ,m+i)- Therefore, 



(27) 



(n,m) 



r„ = LO + W 



Initial conditions (the first of (|24jl ) give 
(28) 



Qo 1 - q 2 0J 2 



P q 2 -to 2 " 
The second and third relations of ()24|) give 



n-l 



(29) | " I = M • I " " I . >( 

Qn / \ Qn-1 



g(l- W 2) 

A/ 



LO 



LO 



-AT 1 ~~ 1 u 



\ q{l-L0 2 ) ~ g (l -lo 2 ) ) 

Let n, and O -1 be the eigenvalues of M. Then one may choose the basis of eigenvectors of A4 
such that 

(30) p n = A ++ n n + A_ + n- n , Q n = A + _n n + A__n~ n . 

The boundary condition \P& = ^A-fc reads in double-indices VPnv-i,m) = ^(o,M-i-m)i 1 - e - 

Pa 



(31) 



N-l 

Qn-i 



M 



N-l 







l-l 






v <3» ; 
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The following par of equations summarizes all the calculations: 



n + 1 = q ;r u l + u- u 1 " q u 



q(l-u 2 ) q(l-u 2 ) 
(32) 

a 2 - O 2 1 - a 2 Q 2 



q (i-n 2 ) q(i-n 2 ) ' 

Here the first equation is the definition of f2 (it is the characteristic polynomial of A4). The 
second equation comes from Q31[) after excluding of Po, Qq via (|28[). The final answer is the 
Bethe Ansatz (two-particles wave function is the superposition of one-particle plane waves) 

(33) * (n>tn) = A ++ Q n cj m + A- + n- n uj m + A+_O n w~ m + A—fl- n uj- m . 
Parameters A± ; ±, related to the eigenvectors of Ad, are defined by 

(34) O iV - 1 A+_ = , = a N - 1 u M - 1 A ++ 
and 

r v A + _ + _ 1 - q 2 o; 2 A_ + + A— _ 1 - q 2 ^ 2 

( j ^+++^-+ ~ q 2 -u; 2 ' + " g 2_^2 • 

The state (|33|) has the evident N <-> M symmetry: 

(36) r a |*) = (w + w -1 )!*) , t p \v) = (n + n- 1 )^) . 

From the point of view of tq, the second relation of (|32j) is the characteristic polynomial, while 
the first one comes from the /3-boundary conditions. 

4.3. p-bosons state. Turn finally to the state with Af = p bosons. Actually, all the results of 
this section were obtained explicitly for p = 3,4, and then conjectured for arbitrary p. 
The wave function f j of the eigenstate 

(37) |*> = Y, vl; /:,/.- , ? .y^y,:'---y./.." 

h<n<-<j p 

is the superposition of the plane waves. In the generic point (no coincidence in ■ ■ ■ ,j P ) the 

wave function is 



(38) .,, ,, = £ K>Ki ■ ■ ■ «3 <S 



-' . • • • 

2 P 



where <r and a' are independent permutations of the set (1, 2, . . . ,p), and n c , m c are related with 
3c by 

(39) j c = Mn c + m c , < m c < M , < n c < N , 
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cf. 1)26)) . Two sets of exponential momenta 
(40) 



uj c = e ik c, n c = e ik '°, c=l,2,...,p 



are solutions of two sets of equations. To write out these equation, we need some extra notations. 
Firstly, let 



(41) 



Ga,b(W}) 



'uj b - quj 

UJ b - UJ a 



Let Ik be a length-A; subsequence of (1, 2, ... , p), I p -k be the compliment subsequence such that 



(42) 

Let then 



IkUlp-k = (l,2,...,p) 



\ 

) 

/ 



(43) Pkiw,co M }) = e n^ M n G bA^}) 

where the sum is taken over all possible subsequences of the length k. E.g. for p = 3 

Pi[{uj,uj m }] = wf f G2iG3i + uj^ 1 G12G32 + G13G23 , 

(44) P 2 u; M }] = ufufG^Gzx + <u,f G21G23 + u%u%GviGvi , 
p [ r, , , ,M -11 _ , ,M, .M. ,M 

The final step: let 

p 

(45) <p(tt|{ W ,w M }) = E(-) fcfiP ~ fc ^[{^ M }] = ^ - JF^Px + ^- 2 P 2 + • • • . 

fc=0 

Then the consistency conditions for the Ansatz 1)38)1 read 

<P(f} a |{^ M }) = 0, 

(46) J V a = l,2,...,p. 

y(u a \{n,n N }) = 0, 

V 

These are the Bethe Ansatz equations for our two dimensional Bose gas. 

The values of ^1,^2,— ,j P m the case when some of j a coincide, as well as the values of (p\) 2 
amplitudes A a ^> , can be defined uniquely. 

The state (|3*H|) provides the eigenvalues of ((11)) 



(47) 

Note, 
(48) 



t |<p) = 1*) (XX 1 ) . r ^i*) = 1*) (fx 1 



\a=l 



\a=l 



17x^2 ' ' ' f^r 



M , M , M 
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(it follows from P p = uof 1 ■ ■ ■ uip 1 ), therefore the mass center of p-particles state moves accordingly 
to 



The one-dimensional limit of (|4fij) must be mentioned. In the case N = 1 our system becomes 
the chain of the length M — 1. Polynomial ^3(u;|{f2, f2}) has a simple structure, all £l a may be 
excluded from (|46|). The resulting equations for u; a are 

what are the usual Bethe Ansatz equations for the quantum chain. 

The last important note is that equations (|46|) provide real momenta k c and k' c of (|40[) if g is 
real. In addition, numerical estimations for not too big lattices and for two particles shows that 
the Ansatz is complete. 

4.4. Sketch derivation of (|46|) . Equations (|46|) may be obtained in the same way as (|32jl . 

One may start with the case when mi < m-2 < • • • < ra p and m p < mi + M. Then the eigenstate 
of T a , r a |*) = |*)r a , r a = ^o;" 1 , is given by 



(50) * mi ,..., mp = E^oC 



m 2 , ,'"p 



The eigenvalue equation for r a provide the way to interpolate $ mi] m 2r .. imp for the larger values 
of m a . For instance, when m p ~ mi + M, the eigenvalue equation modifies as 

(51) r a * mi) ..., mp = g' 5m f' mi+M - 1 ^ mi -i,..., mp + ^ mp ' mi+M ^ mi ,...,m p -i + all the rest . 

Extra g-factors produce some linear transformation of the amplitudes A a of l I / mii ... imj ,| mp > mi +M 
in the some way as in eq. (|29|). Repeating this procedure for m p ~ m2 + M, m p ~ m3 + M etc., 
one comes to 



(52) * mi , m2 ,..., mp+ M = ^4 1} < a 



<T2 °"p ' 



where is a linear combination of A^: 



(0) 



(53) AW = £ 4° 

cr' 

cf. 1|29|) . The Ansatz (|38|) corresponds to (|52|) in the basis of eigenvectors of A 7 !. The miracle of 
the exact integrability is that p\ x p\ matrix A4 has only p eigenvalues: 

(54) det(n-M) = WiflHuj,^ 1 })^-^ 

The second equation in (|46|) . providing the N M symmetry of the Bethe Ansatz equations, 
solves the boundary condition analogously to (pTT|) . 
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s. sergeev 
5. Conclusion 



This paper presents a model of two-dimensional Bose gas and the conjecture for the equations 
(|46Jl describing its eigenstates. The model has the main features of a physical model. The states 
are described by the system of plane waves - at least the notion of the momenta is well defined. 
In particular, the sum of cosines of momenta is a candidate for the Hamiltonian. Equations 
(|46|). in the case when the number of bosons p is big, are the subject of further investigations. 
We believe, this model may have some interest for the theory of integrable systems and for the 
condensed matter physics. 

Acknowledgements. I would like to thank M. Batchelor, V. Bazhanov, M. Bortz, X.-W. Guan 
and V. Mangazeev for valuable discussions. 
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